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NOTATION 


One  half  the  average  molecular  speed 
Velocity  vector  of  molecules  relative  to  U 
Specific  heat  at  constant  pressure 
Specific  heat  at  constant  volume 
Internal  energy  per  unit  mass 
Body  force  per  unit  mass  potential 
Idem  factor,  I = ii  + jj  + kk 

Outward  drawn  unit  vector  normal  to  the  surface 
element  dS 

Static  pressure 

Energy  flux  vector 

Heat  addition  rate  per  unit  mass 

Universal  gas  constant  divided  by  molecular  weight 

Closed  surface  bounding  V 

Square  root  of  the  quantity  in  parenthesis 

Vector  tangent  to  the  surface  element  dS 

Absolute  static  temperature 

Flow  velocity  vector 

Arbitrary  volume 

Flow  property  value  at  the  surface  element  asso- 
ciated with  those  molecules  going  out  of  V 

Flow  property  value  at  the  surface  element 
associated  with  those  molecules  going  into  V 

Average  distance  of  molecular  collisions  from  the 
surface  element,  two-thirds  of  the  mean  free  path 

Viscosity 

3.14159 


iv 


p 


Density,  mass  per  unit  volume 


t Average  transit  time  of  molecules  from  collision 

point  to  the  surface  element,  t = A/a 


Special  Symbols 

3(  ) Partial  time  derivative  of  quantity  in  parenthesis 


( ) ' Property  of  molecules  moving  out  of  V through 

the  surface  element 

( )"  Property  of  molecules  moving  into  V through 

the  surface  element 

( )c  Conjugate  of  ( ) 


ABSTRACT 


Starting  with  the  frictionless  flow 
i equations,  governing  equations  for  viscous 

compressible  flow  are  developed  from  a phys- 
ical molecular  transport  model  by  forming  a 
first  perturbation  of  the  average  flow  prop- 
erty value  transported  across  a surface  ele- 
ment by  two  distinct  groups  of  molecules. 
This  results  in  a continuity  equation  with 
viscous  correction  terms. 


ADMINISTRATIVE  INFORMATION 
This  investigation  was  authorized  and  funded  by  the 
Naval  Air  Systems  Command  (AIR-320D)  under  Project  Element 
62241N,  WF  41.421.091,  Work  Unit  1-1600-078,  and  was  accom- 
plished during  the  months  of  April,  May,  and  June  1978. 


INTRODUCTION 

The  frictionless  flow  equations  may  be  written  in  terms 
of  a volume  integral  equal  to  a surface  integral,  where  the 
surface  integral  integrand  accounts  for  flow  of  some  physical 
property  across  a surface  element.  The  molecules  crossing  the 
surface  element  fall  into  two  groups,  those  with  a component 
of  relative  velocity  in  the  same  direction  as  the  surface 
element  unit  vector  n and  those  with  a component  opposite  to 
n.  By  applying  the  simple  physical  molecular  model  in  which 
some  average  flow  property  is  transported  from  the  last  aver- 
age collision  point  across  the  surface  element  to  these  two 
distinct  groups  of  molecules,  correction  terms  are  obtained 
for  the  frictionless  flow  continuity  equation  as  well  as 
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correction  terms  for  the  momentum  and  energy  equations. 

The  correction  terms  for  the  continuity  equation  are  found 
to  be  essential  for  the  solution  of  certain  physical 
problems . 


DEVELOPMENT 

The  frictionless  flow  equations  may  be  written  using 
dyadic  and  vector  notation  in  terms  of  volume  and  surface 
integrals-*-* 


f 3(p)dV  = - / n • ( pU)dS 


/ ta(pU)-pVF]dV  = - / n • ( pUU+pnn)dS 


/ f 8 ( pe+pU* U/2-pF ) -pq] dV 


= - ; n.  ( pe  + pU»U/2-pF+p)UdS 


where  9()  denotes  the  partial  derivative  with  respect  to 
time,  VF  is  the  body  force  per  unit  mass  due  to  the  potential 
F,  and  -F  is  the  potential  energy  per  unit  mass.  These  equa- 
tions neglect  molecular  transport  effects  so  there  is  no 
viscosity  or  thermal  conductivity.  Thus,  in  a volume  element 
containing  the  surface  element  ndS,  the  molecular  velocity 
distribution  is  locally  Maxwellian.  Then,  half  of  the 


*A  complete  listing  of  references  is  given  on  page  17, 
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molecules  with  density  p'  will  have  velocity  c'  relative 
to  u so  that 

n.c'  > 0 (4) 

and  the  remaining  molecules  with  density  p"  will  have 

n.c"  <0  (5) 

Then  for  Equation  (1) 

pU  = p' (U+an)  + p" ( U-an ) (6) 

identifies  the  two  groups  of  molecules  crossing  ndS,  P'  from 
inside  V going  out  and  P”  from  outside  V going  in,  where 

P'  = P"  = P/2  (7) 

and 

a = SQR(2RT/tt)  (8) 

In  Equation  (8),  a is  found  to  be  one-half  the  average 
molecular  speed  from  kinetic  gas  theory,^  R is  the  universal 
gas  constant  divided  by  the  molecular  weight,  and  T is  the 
absolute  static  temperature. 

Using  p = pRT  in  Equation  (2) 

(pUU+pnn)  = ( pUU+pRTnn) 

= [p1 (U+an) (U+an) +p"( U-an) (U-an) ] (9a) 

identifies  the  two  groups  of  molecules  crossing  ndS  and  the 
velocity  associated  with  each  group,  where  it  is  understood 
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that  paa  must  be  multiplied  by  tt/2  when  it  occurs  to 
account  for  root  mean  square  averaging  and  obtain  the  pRT 
term.  Expanding  Equation  (9a)  and  using  the  idem  factor^  I 
to  replace  nn,  because  n»nn  = n»I  = n,  one  finds 

( pUU+pnn ) = p ' (UU+anU+aUn+RTI) 

+ p" (UU-anU-aUn+RTI)  (9b) 

Using  p = pRT , e = cvT,  and  cp  = R + cv  in  Equation  (3) 

p(cpT+U*U/2-F)U  = p’ (cpT+U«U/2-F) (U+an) 

+ p" (cpT+U*U/2-F) ( U-an ) (10) 

identifies  the  two  groups  of  molecules  crossing  ndS  and  the 
energy  associated  with  each  group. 

To  form  the  perturbation,  assume  some  average  flow  prop- 
erty value,  e.g.,  X',  is  acquired  by  the  p'  molecules  at  their 
last  average  collision  point  before  crossing  ndS. 

Then,  if  X is  the  average  flow  property  at  the  surface  element 

X'  = X - ( i/a) (U+an). V(X)  (11a) 

is  the  average  flow  property  value  acquired  by  the  p'  mole- 
cules and  transported  across  the  surface  element,  where  X 
is  the  average  distance  of  the  collision  point  from  dS  when 
n«U  = 0.  Also,  X is  two-thirds  of  the  mean  free  path.  The 
ratio  (X/a)  is  an  average  characteristic  time,  e.g.,  t. 
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between  molecular  collisions  and  passage  of  the  molecules 
through  dS.  Thus,  Equation  (11a)  may  also  be  written 


X'  = X - x( U+an ) • VX  (lib) 

Similarly,  if  Y"  is  an  average  flow  property  value  of  the  p" 
molecules 

Y"  = Y + t ( an-U) • VY  (12) 

The  time  averaging  of  molecular  motion  in  a volume  ele- 
ment provides  the  basis  for  the  perturbation  model  and  this 
process  should  be  examined  briefly.  At  standard  conditions 
there  are  approximately  2.7  x 10-*-9  molecules  in  a cubic  centi- 
meter of  gas.  The  dimension  of  X is  approximately  4.3  x 10-® 
cm  and  a is  approximately  2.2  x 104  cm/s,  hence,  t = 2 x 10-1® 
seconds.  Then,  as  the  molecules  enter  and  leave  a cubic  vol- 
ume element  of  dimension  2A,  there  will  be  approximately 
6 x 8.6  x 8.6  x 2.2  x 2.7  x 10^  = 2.6  x 1074  discrete  values 
of  a particular  frictionless  flow  property  (p,U,p,etc.)  per 
second.  If  the  macroscopic  time  element  associated  with  the 
flow  properties  in  the  frictionless  flow  equations  is  taken 
as  one  microsecond,  then  during  that  time  element  there  will 
be  approximately  2.6  x 10®  discrete  values  of  the  frictionless 
flow  property  that  form  the  average  value  associated  with 
the  location  of  the  cubic  volume  element  at  that  particular 
time.  During  this  microsecond,  approximately  2.2  x 107  mole- 
cules will  enter  a side  of  the  cubic  volume  element  and 


penetrate  to  its  center;  thus  transporting  their  flow  property 
value  across  a surface  element  at  the  center  of  the  volume 


element.  The  flow  property  value  transported  across  the  sur- 
face element  by  these  penetrating  molecules  may  be  different 
from  the  average  frictionless  flow  property  value  of  the  volume 
element,  which  is  macroscopical ly  flowing  through  the  surface 
element.  The  perturbation  model  attempts  to  account  for  this 
difference.  To  be  consistent,  the  perturbation  model  will 
be  applied  to  the  surface  element  flux  of  all  three  friction- 
less flow  equations. 

Substituting  Equations  (6),  (lib),  and  (12)  in  Equation 
(1)  forms  the  continuity  perturbation, 

/ 9( p)dV  = - / n* [ p' ( U+an)+p" (U-an) 

-x( U+an)*  7( p' U+p'an)+x(an-U). 7( p"U-p"an) ]dS 
= - / n • [ pU-xU*  7( pO) -xan»  7( pan ) ] dS  (13a) 

where  X = p'(D+an)  and  Y = p"(U-an)  has  been  used. 

Expanding  the  last  term  of  Equation  (13a)  one  obtains  a com- 
ponent n«(n.7n).  The  derivative  n«7n  will  be  a vector  t 
which  is  tangent  to  dS.  Thus,  n*(n*7n)  = n*t  = 0.  With  this 
and  n*n  = 1,  Equation  (13a)  may  be  rewritten 

f 9(p)dV  = - / n • [ pU- xU • V ( pU ) - A 7 ( pa ) ] d S (13b) 


The  last  two  terms  of  Equation  (13b),  dotted  by  n,  represent 
a decrease  of  the  frictionless  flow  mass  flux  n»(pU)  crossing 
dS  due  both  to  a shift  of  the  last  two  average  collision  points 
by  the  frictionless  flow  velocity  and  to  a molecular  mass 
flux  gradient.  In  the  last  term  of  Equation  (13b),  the  X may 
be  replaced  by  viscosity,  y = paX.  With  this  and  the  diver- 
gence theorem.  Equation  (13b)  becomes 

f {3( p)+V.( pU)-V*[xU.V( pU)+yVln( pa)] }dV  = 0 (14) 

Because  the  volume  is  arbitrary,  the  integrand  of  Equation 
(14)  must  be  zero.  Thus,  the  continuity  equation  becomes 

3 ( p)  + V.(pU)  = V* ( TU.V( pU)+yVln( pa) ] (15) 

where  the  terms  on  the  right  provide  corrections  to  the 
macroscopic  frictionless  flow  continuity  equation  due  both  to 
a shift  of  the  last  two  average  collision  surfaces,  caused  by 
the  flow  velocity,  and  to  the  molecular  mass  transport 
gradient . 

Substitute  Equations  (9a),  (9b),  (lib)  and  (12)  in  Equa- 
tion (2)  to  form  the  momentum  perturbation 
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/ [ 3(  pU)-p7F]dV 


= - f n*{p* (U+an) ( U+an  ) + pB  (U-an) (U-an) 

- t( U+an ) • V [ P' (U+an ) ( U+an ) ) 

+ x( an-U) • V ( p" (U-an) (U-an) ] )dS 

- - f n. [ pUU+pI-TU*V( puD+pI ) 

-An*7( panU+paUn) ]dS  (16) 

Using  the  divergence  theorem 

/ [9( pU)+tf*( pUU)+Vp-pVF-V. [rU.V(pUU) ] 

-7(xU»7p)]dV  = / n • [ An • V( panU+paUn ) ] dS  (17) 

The  surface  integral  in  Equation  (17)  contains  the  symmetric 
dyadic  (panU+PaUn)  whose  value  depends  on  the  orientation  of 
n.  Carrying  out  the  differentiation  for  n*7,  the  surface 
integral  integrand  becomes 

n • [ An  »7( panU+paUn ) ] 

= n.A{(n[n.7(paU)]+t(paU) 

+ [n «V( paU) ] n+( paU) t ) } (18) 


The  dyadic  on  the  right  side  of  Equation  (18)  is  also 
symmetric.  Using  n«n  = 1 and  n*t  = 0,  this  may  be  rewritten 

n • [ An  • V ( panu+  paUn ) ] 

= n • A [ V( paU ) + [n«  V( paU)  n+( paU) t ] (19) 

The  dyadic  on  the  right  side  of  Equation  (19)  must  be  sym- 
metric to  preserve  equilibrium  of  the  volume  element,-*  thus 
the  last  terms  must  be  replaced  by  [V(paU)]c,  which  is  the 
conjugate  of  V(paU).  The  normal  component  n»fn«V(paU)]  is 
found  to  be  equal  to  the  normal  component  n • ( n • [7( paU) ] c ) by 
merely  expanding;  thus  it  is  the  tangential  component  of 
n»(V(paU))c  which  must  be  equal  to  n*(paU)t  to  preserve  the 
equilibrium  of  the  volume  element.  Other  authors5  showed 
the  requirement  for  a symmetric  dyadic  using  the  principle 
of  angular  momentum  conservation,  so  that  this  view  could 
be  used  instead  of  preserving  equilibrium  of  the  volume 
element.  Perhaps  the  opposite  point  of  view  would  be  to 
assume  that  a volume  element  is  not  in  equilibrium  and  has 
started  to  spin  up.  Then,  in  a microsecond,  the  molecules 
within  the  volume  element  would  flux  back  and  forth  about  a 
thousand  times,  transporting  tangential  momentum  which  opposes 
the  spin  to  reestablish  the  volume  element's  equilibrium. 
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In  any  case,  the  dyadic  must  be  symmetric.  With  Equation  (19) 
and  the  divergence  theorem.  Equation  (17)  becomes 

/ { 9 ( pU)  +V • ( pUU) +Vp-pVF-V • [ tU • V( pUU)] 

-V( tU« Vp) - V» ( AV( paU) +X [ V( paU) ] c ) }dv  = 0 (20) 

Because  the  volume  is  arbitrary,  the  integrand  of  Equation 
(20)  must  be  zero.  Using  y,  the  momentum  equation  becomes 

9 ( pU)  + V*(PUU)  + 7p  - pVF  = V • [ xU*  V( pUU) ] 

+ V(tU.7p)  + [(p/pa){V(paU)  + [V(paU)]c}]  (21) 

where  the  first  two  terms  on  the  right,  involving  t,  provide 
a correction  due  to  a shift  of  the  last  average  collision 
surfaces  by  the  flow  velocity.  The  last  terms  provide  a 
correction  due  to  the  molecular  momentum  transport  gradient. 
The  dyadic  V( paU)  and  its  conjugate  may  be  expanded  to  show 
the  correction  that  is  due  strictly  to  viscosity  and  velocity 
gradient 

( y/pa ) (7( paU)  + [7( paU) ] c } 

= y[7U+(7U)c]  + y ( ( Vln  pa ) U+U( 7lnpa ) ] (22) 
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Using  Equations  (10),  (lib),  and  (12)  in  Equation  (3) 


to  form  the  energy  perturbation 

f { 3 ( pe+pU »U/2-pF ) +V • [ ( pe+pU« U/2-pF  + p)U]-pq}dV 
= f n •{ t ( U+an ) • V [ p ' ( cpT+U»U/2-F) ( U+an ) ) 
+T(U-an).V[p" (cpT+U*U/2-F) (U-an) ] }dS 
= / n.{xU*V[p(cpT+U»U/2-F)U] 

+X7[pa(cpT+U*U/2-F) ] }dS  (23) 

With  the  divergence  theorem,  Equation  (23)  may  be  written 
as  a volume  integral  equal  to  zero,  and,  since  the  volume  is 
arbitrary,  the  integrand  is  zero.  Then,  using  X = y/(pa), 
the  energy  equation  may  be  written 

3 (p [cvT+U*U/2-F]  ) + V • [p ( CpT+U« U/ 2-F) U]  - pq 

= ^»{tU»V[p( CpT+U* 0/ 2-F ) U] +u V( CpT+U* U/2  -F) 

+ U ( CpT+U*U/2-F ) V ( In ( p a ) ] } (24) 

DISCUSSION 

The  significant  result  obtained  from  this  development 
is  the  continuity  equation  with  molecular  transport  correc- 
tion terms,  Equation  (15).  These  terms  correct  the  macro- 
scopic frictionless  flow  continuity  equation  and  are  the  same 


order  of  magnitude  as  the  correction  terms  for  the  momentum 


and  energy  equations.  It  should  also  be  noted  that  the 
momentum  and  energy  equations,  given  by  Equations  (21)  and 

I 

(24),  provide  new  correction  terms  due  to  a shift  of  the 
molecular  collision  surfaces,  on  both  sides  of  the  surface 
element  dS,  by  the  flow  velocity.  When  viewed  in  the  surface 
integrals,  all  of  these  correction  terms  may  be  physically 
interpreted  as  perturbations  to  the  locally  Maxwellian  mole- 
cular velocity  distribution  function.  The  perturbations  re- 
sult from  the  collision  surfaces  being  situated  on  both  sides 
of  the  surface  element  dS,  rather  than  coincident  with  it. 

Two  simple  examples  will  illustrate  the  nature  of  these 
correction  terms.  First,  consider  a thermally  insulated 
container  of  gas  at  rest,  with  uniform  pressure,  without  body 
forces,  and  with  a thermally  insulated  partition  so  that  the 
two  compartments  of  the  container  are  at  different  tempera- 

t 

tures.  Now,  if  the  partition  is  removed,  there  will  be  un- 
steady effects  due  to  molecular  motion  which  initially  trans- 
port mass  and  energy  across  the  surface  elements  where  the 
partition  was  located.  For  the  two  compartments,  let  T2  < T^, 
and  because  the  pressures  are  equal  p^  < P2»  The  net  mass 
flux  from  1 to  2 will  be 

( Piai~P2a2 

* [ P/SQR(  2Rtt  ) ] [1/SQR(T1)-1/SQR(T2)]  (25) 

it 

i 
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The  net  energy  flux  from  1 to  2 will  be 


( p1a1cpT1-P2a2cpT2)/2 

= [pcp/SQR(2RTT)]  [SQR(T1)-SQR(T2)J  (26) 

Thus,  near  the  surface  element  there  will  be  a time  rate  of 
change  of  density  8(  p)  even  though  the  macroscopic  flow 
velocity  is  zero.  This  effect  is  included  in  Equation  (15). 
Since  these  two  net  effects  occur  at  the  molecular  level 
within  a few  mean  free  paths,  a net  pressure  gradient  will 
be  established  across  the  surface  element  which  starts  to 
drive  the  momentum  equation.  Of  course  the  system  without 
the  partition  will  eventually  reach  a new  state  of  thermo- 
dynamic equilibrium. 

As  a second  example,  consider  the  flow  near  the  center 
of  a long  thermally  insulated  pipe  filled  with  gas,  without 
body  forces,  which  is  conducting  heat  in  a steady  state  con- 
dition. Taking  Equation  (15)  back  to  a surface  integral,  the 
macroscopic  mass  flux  will  be 

pU  = xU  *V( pU ) + pVln(pa)  (27) 

which  is  exactly  balanced  by  the  small  microscopic  mass  flux 
in  the  -U  direction. 
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Taking  Equation  (24)  back  to  a surface  integral,  the  energy 
flux  vector  Q will  be 

Q = p(CpT+U* U/2 )U  - TU«7[ P( CpT+U* U/2 ) U] 

- p7(cpT+U*U/2)  - U(cpT+U*U/2) Vln( pa)  (28) 

This  equation  shows  that  the  energy  is  transported  macro- 
scopically  by  pU  as  well  as  microscopically  by  the  molecular 
activity.  Using  Equation  (27),  the  energy  flux  vector  may 
be  written 

Q = - [ tU • 7 ( pU) +yVln ( pa ) ] [ tU*  7(  cpT+U»  U/2 ) ] 

- y7 ( CpT+U« U/2 ) (29) 

The  principal  term  contributing  to  the  energy  flux  vector  in 
Equation  (29)  is  -uV(CpT).  If  the  body  force  potential  F 
were  not  zero,  the  energy  flux  vector  Q would  include  this 
effect  from  Equation  (24).  Thus,  if  a long  pipe  were  used 
to  experimentally  determine  the  coefficient  of  thermal  con- 
ductivity, the  result  would  be  affected  slightly  by  the 
orientation  of  the  pipe  to  the  potential  gradient.  In  this 
case,  the  approximate  net  energy  flux  would  be 

n*Q  * -yn»V(CpT-F)  (30) 

The  correction  terms  involving  tU  in  Equations  (15), 
(21),  and  (24)  will  become  significant  in  high-speed  flows. 
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such  as  found  in  ballistic  reentry  and  shock  wave 
calculations.  This  occurs  when  U approaches  the  same  order 
of  magnitude  or  is  greater  than  a,  so  that  tu  is  several  or- 
ders of  magnitude  larger  than  in  the  two  examples  given  above. 

CONCLUSIONS 

The  governing  equations  for  viscous  compressible  flow 
have  been  derived  by  using  a simple  molecular  transport  model 
to  form  an  approximate  perturbation  of  the  frictionless  flow 
equation.  The  molecular  transport  correction  terms  in  the 
continuity  equation  are  of  the  same  order  of  magnitude  as  the 
correction  terms  in  the  momentum  and  energy  equations.  It 
appears  that  all  of  these  correction  terms  should  be  con- 
sidered initially  for  the  proper  analysis  of  diverse  flow 
processes . 
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